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Abstract 

We apply a discrete version of the methodology in [12] to obtain a recursive asymp- 
totic expansion for E[/i(H^)] in terms of Poisson expectations, where is a sum of 
independent integer- valued random variables and his a polynomially growing function. 
We also discuss the remainder estimations. 
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1 Introduction and main result 

It should be noted in the first place that the notation used in this paper is similar as in 
|12j . however, their meanings are different since we here consider discrete random variables. 
Stein's method for Poisson appoximation has been introduced by Chen [8] . Let Z be an N- 
valued random variable (N-r.v.), then Z follows the Poission distribution with parameter 
A if and only if the equality E[Zf{Z)] = XE[f{Z + 1)] holds for any function / : N ^ M 
such that both sides of the equality are well defined. Based on this observation, Chen has 
proposed the following discrete Stein's equation: 

(1) xf{x)-Xf{x + l)=h{x)-Vx{h), x€N 

where 'P\{h) is the expectation of h with respect to the A-Poisson distribution. If X is an 
N-r.v., one has E[/i(X)] - Vx{h) = ¥.[Xfh{X) - fh{X + 1)] where fh is a solution of ((!]) 
and is given as 



( 1 M \ ^ 

(2) M^)=^-^E^i''(^-'P>^(^))- 



l=X 
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The value //i(0) can be arbitrary and is not used in calculations in general. 

Stein's method has been adopted for Poisson approximation problems since ^ in a 
series of papers such as |T], [5], [3j among many others, one can also consult the monograph 
[6] and the survey paper [10]. In particular, Barbour [3] has developed, in parallel with the 
normal case p], asymptotic expansions for sum of independent N-r.v.s and for polynomially 
growing functions. The asymptotic expansion problem has also been studied by using other 
methods such as Lindeberg method (e.g. 0). 

In this paper, we address this problem by the zero bias transformation approach. 
Similar as in Goldstein and Reinert [H], we introduce a discrete analogue of zero bias 
transformation (see also [9]). Let X be an N-r.v. with expectation A. We say that an 
N-r.v. X* has Poisson X-zero biased distribution if the equality 

(3) IE[^/(^)] = AE[/(X* + 1)] 

holds for any function / : N ^ M such that the left side of ^ is well defined. The distri- 
bution of X* is unique: one has P(X* = x) = {x + 1)P(X = x + 1) / X. Combining Stein's 
equation ([T|) and zero bias transformation ([3]) , the error of the Poisson approximation can 
be written as 

(4) nKX)] - Vx{h) = mux* + 1) - h{X + 1)]. 

A first order correction term for the Poisson approximation has been proposed in ^ by 
using the Poisson zero bias transformation. 

Recall the difference operator A defined as A/(x) = f{x + 1) — f{x). For any x G 
N* := N \ {0} and any n € N, one has A(^) = {^^- If / and g are two functions on N, 
then 

A(/(x)5(x)) = /(x + l)A5(x) + 5(x)A/(x). 

We have the Newton's expansion (^ ThmS.l]), which can be viewed as an analogue of 
the Taylor's expansion in the discrete case. For all x, y € N and € N, 

^ / \ 

/(x + y) = ^rjAV(x)+ Yl A^+V(x + ii). 



Let us introduce the following quantity, where we use the same notation as in [T2J , but 
its meaning is changed. For any N-r.v. Y and any /c S N such that E[|y['^] < +oo, denote 

by 



(5) 



m 



Y 



E 



kl [Y]k 



where [Y]k is the /c*^ factorial moment of Y. Let X and Y be two independent N-r.v.s 
and / : N ^ M such that A''f{X) and A^/(X + Y) are both integrable, then 



N 



(6) 



E[/(X + Y)]=Y, m^EiA'^fiX)] + 5^(/, X, Y). 



k=0 
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where 

(7) SN{f,X,Y)=E[ A^^'f{X + h] 

0<ji<-<jN+i<Y 

We introduce the discrete reverse Taylor formula. Once again, the following result is 
very similar with [121 Prol.l], however, with different significations of notation. 

Proposition 1.1 (discrete reverse Taylor formula) With the above notation, we have 

(8) E[f{X)]=Y,i-'^)' E rn?E[A\'\fiX + Y)]+er,if,X,Y) 

d>0 3emd:,\3\<N 

where 

(9) ^^(/,x,y):=-E(-l)'^ Yl rn^Y^^N-\J\iA\'\f,X,Y), 

d>0 JGN^,|J|<Ar 

for any integer d > 1 and any J = (ii)f=i ^ |J| = ji + ■ ■ ■ + Jd cii^'d rriy"^'-* : = 
ruy^^ ■ ■ ■ mY'^\ and by convention, = {0} with |0| = 0, ruy^ = 1. 

Consider now a family of independent N-r.v.s Xi{i = 1, ■ • - n) with expectations Aj, 
which are "sufficiently good" in a sense we shall precise later. Let W = Xi + • • • + X„ 
and denote Xw ■= ^W] = Ai + • • • + A„. Let W^'^ = W - Xi and X* be an N-r.v., 
independent of and which has the Poisson Xj-zero biased distribution. Finally, let I 
be a random index valued in {1, • • ■ , n} which is independent of [Xi, • • • , X*, • • • , X*) 
and such that P(/ = i) = Aj/A^y for any i. Then, similar as in yilj, the random variable 
W* := + X} follows the Poisson W -zero biased distribution. 

We give below the asymptotic expansion formula in the Poisson case. 

Theorem 1.2 Let € N and p > 0. Let /i : N — > M 6e a function which is of 0{xP) at 
infinity and Xi {i = 1, ■ ■ ■ , n) be a family of independent N-r.v.s having up to {N +p+ 1)*'^ 
order moments. Let W = Xi + - ■ ■ +Xn and Xw = E[W]. Then K[h(W)] can be written as 
the sum of two terms Cj^ih) and e^ih) such that Co{h) = Vx^ih) and eo{h) = lE,[h{W)] — 
Vx^ih), and recursively for any N >1, 

(10) CM{h) = C,{h) + j2X.Y.^-lY~'T. m(?;)(m(?!)-m^^;))C^„|j|(AlJlA(x + l)), 

i=l d>l JeN^,|J|<7V 

eM{h)=jZx{Y.i-^f-' E -rKf -m?;V^-lJ|(Al^lA(x + l)) 

«=1 d>l JeN^,|J|<Af 
^ ' N 

+ V mf. EN-kiA'^fhix + l),W^'\Xi) + 6N{fh{x + 1), W^'\X*)\ , 

k=0 

where for any integer d > 1 and any J £ N^, G N,, denotes the last coordinate of 3, and 
J° denotes the element in N!f~^ obtained from 3 by omitting the last coordinate. 
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Remark 1.3 In view of the similarity between the above theorem and pL2, Thml.2], 
which has also been shown by the two papers [21 [3] of Barbour, the following question 
arises naturally: can we generalize the result to any infinitely divisible distribution? 

2 Several preliminary results 

In this section, we are interested in some properties concerning the function h and the 
associated function //j. Compared to the normal case, we no longer need differentiability 
conditions on h in Theorem 11.21 and shall concentrate on its increasing speed at infinity. 
This makes the study much simpler. 

We begin by considering the modified Stein's equation on N*: 
(12) xf{x) - Xf{x + 1) = h{x), X G N,. 

The above equation may have many solutions, one of which is given by 

7f \ (x - 1)! ^ X\^.^ 

i=x 

A general solution of (jl2p can be written as fh{x) + C(a; — 1)!/A^, where C is an arbitrary 
constant. However, when h is of polynomial increasing speed at infinity, fh is the only 
solution of (jl2p which has polynomial increasing speed at infinity. 

In order that the function //^ is well defined, we need some condition on h. Denote by 
£x the space of functions h on N* such that, for any polynomial P, we have 

i>l 

Clearly £\ is a linear space. We list below some properties of £\. 
Proposition 2.1 The following assertions hold: 

1) for any Q G M[x,x^^] and any h ^ £x where denotes the set of Laurent 
polynomials on M, we have Qh G £x; 

2) for any h G £x, Ah G £x and fh G £x; 
Proof. 1) is obvious by definition. 

2) Let hi be the function defined as /ii(x) := h{x + 1). If P is a polynomial, then 

i>l ■ j>2 j>2 
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since h (z £\. Therefore, Ah € £x. We next prove the second assertion. For any arbitrary 
polynomial P, there exists another polynomial Q such that, for any integer i > 1, Q{i) > 
Yl)=i ^^j^. Therefore 

a>l i=a o>l i=a 

= E (E^)^iM^)i < E^w^iM^)! < 

i>l a=l ' i>l 

which implies that ^ £x. □ 

For any function /i G we define t(/i) : N ^ M such that 

T{h){x) =h{x + l)/{x + l). 

Note that for any integer A; > 1, one has T^{h){x) = xlh{x + k)/{x + k)\. The proof of 
Proposition 12.11 shows that r is actually an endomorphism of £x. 

Lemma 2.2 Let h ^ £\. Then 

(13) U^h^{x) = Jh{x + l)/x. 

Proof. Let u[x) = fh{x + ^)/x. Dividing both sides of (jl2p by x and then replacing x by 
X + 1, we obtain fh{x + 1) — ^fh{x + 2)/{x + 1) = T{h){x), or equivalently, 

(14) xu{x) - \u{x + 1) = T{h){x). 

Since fr{h) is the only solution of (fHI) in £'a, the lemma is proved. □ 

Corollary 2.3 Let peR. If h{x) = 0{xP), then fh{x) = 0(xP-^). 
Proof. First of all. 



- ^ il ^ i\/x\ ~ ^ ii- x)\ 

1>X 1>X 1>X 



e 



Therefore, when p <0, one has 



rl °° A* 
xf,{x) = -Y,-^Hi)=0{xn 

i=x 

since < if x < i. Hence fh{x) = 0(x^~^). The general case follows by induction on 
p by using □ 
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We now introduce the function space: for any p € M, denote by Tip the space of all 
functions /i : N ^ M such that h{x) = 0{x^) when x — > oo. In the following are some 
simple properties of Tip, their proofs are direct. 

Proposition 2.4 1) For any p > and any h G Tip, the restriction of h on lies in 
riA>o^A- 

2) If h G Tip, then also are h{x + 1) and Ah. 

3) If h G Tip and g E Tig, then gh E Tip+q. 

The following proposition is essential for applying the recursive estimation procedure. 
Proposition 2.5 Let p>0. If h e Tip, then fh € Tip. 

Proof. Note that fh coincides with on N=k where h = h — T^x(h). Since h £ Tip, also is 
h. Then Corollary O implies fh{x) = 0{xP-^) = 0{xP). □ 



3 Proof of the main result 

In this section, we give the proof of Proposition 11.11 and of Theorem [L2l which are essen- 
tially the same with the ones of [12, Propl.l, Thml.2] in a discrete setting. 

Proof of Proposition 11.11 We replace E[AI-'I/(X + Y)] on the right side of <^ by 

N-\J\ 

and observe that the sum of terms containing 5 vanishes with eN{f,X,Y). Hence the 
right side of ([8]) equals 

Af-|J| 

^(-1)^ ^ m? Y: mPE[A\^\^^fiX)] 

d>0 JeN^,|J|<Ar k=0 

If we split the terms for k = and for 1 < k < N — \J\ respectively, the above formula 
can be written as 

d.>0 JeN^,|J|<Ar d>0 JeN^,|J|<Ar k=l 

We make the index changes J' = (J, k) and n = d + 1 in the second part of ()15p and find 
that it is nothing but 

u>l J'eN", |J'|<Af 
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By taking the sum, it only remains the term of index d = in the first part of (jl5p . which 
is equal to E[/(X)]. So the lemma is proved. □ 



Proof of Theorem 11.21 We prove the theorem by induction on A^. The case where = 
is trivial. Assume that the assertion holds for 0, • • • , — 1. 

Since h G TCp, by Lemma [231 and Proposition[23]2), for any k £ {!,■■■ , N}, A''fh{x + 1) e 
Hp-i C Hp. Therefore C^^ki^^ fh{x + 1)) and e7v-A:(A^/ft(x + 1)) are well defined and 

n^'^hiW + 1)] = CN-ki^^'Ux + 1)) + e^_fc(AV/»(x + 1)). 

We now prove the equality E[/i(VF)] = CN{h) + eAr(/i). Recall that for any i G {1, • • • , n}, 
X* follows the Poisson Xj-zero biased distribution and is independent of W'^'^^ = W — Xi, 
I is an independent random index such that P(I = i) = Xi/Xw, and W* = W^^^ + Xj. So 
E[h{W)] - Co{h) is equal to 

n 

XwnfhiW* + 1) - MW + 1)] = Y1 A, (e[A(VF« + X* + 1)] - E[fh{W + 1)] 

i=l 

where, by using ([6]), 

N 

If 

k=0 



E[A(VF» + X* + 1)] = Y.m^m'hiW^'^ + 1)] + SNifhix + 1), X*). 



By replacing E[A'^//i(I^(*) + 1)] in the above formula by its {N — fc)**^ order reverse Taylor 
expansion, we obtain that E,[fh{W^^^ + X* + 1)] equals 



N 

k=0 ' '-d>0 J6N^ 
|J|<Af-fc 



+SNifhix+l),W^'\X*). 



Note that the term of indices A; = d = in the sum is E[fh(W + l)]. Therefore, E[fh{W'^'^ + 
■^i + 1)] ~ + 1)] is the sum of the following three terms 

N 

(16) E-xlE(-l)' E rn^Sm^'^^'MW + l)], 

k=l d>Q JeN^, |J|<Af-fc 

(17) Y(-^^' E "^?>[ai-^ia(w^+i)], 

d>i jmi,\j\<N 

N 

(18) J^m5J]e^_fc(AV/^(x + l),W^'\Xi) + 5^(A(x + 1),W^'\X*). 

k=0 

By interchanging summations and then making the index changes K = (J, /c) and u = d+1, 
we obtain 

(USD = E(-l)""' E ^T^rnfpElA^'^^MW + l)]. 

u>l KeN^,|K|<Af 
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As the equality = m^. holds for any J, (fT6l) + (fT7Il simplifies as 

d>l JGN^,|J|<Af 

By the hypothesis of induction, we have 

E[A\'\MW + 1)] = C^_|j|(AlJ| A(x + 1)) + ejv_|j|(Al-^lA(x + 1)), 

so the equahty E[/i(Vl^)] = C^ih) + ejyih) holds with C^ih) and e^ih) being defined in 
(UnD and (HH). 

□ 



4 Error estimations 



In this section, we concentrate on the remainder e]\f{h) in the asymptotic expansion. The 
following quantity will be useful. Let p >0. For h G Tip and € N, we define 

(19) - ' 



N,p ■■= sup 



xP 



which is finite by Proposition 12.41 2) 



Lemma 4.1 Let N ^ N, k e {0, ■ ■ ■ , N} and p > 0. Let X he an N-r.i). with p*'' order 
moment, Y be an N-r.w. independent of X and having {N — A; + 1 + pY^ order moment. 
Then, for any f G Tip, the following inequalities hold: 



(20) |5jv-fc(AV(2; + l),X,y)| <max(2P-\l)||/|U,p(E[X^']m 



{N-k+l) (N~k+l),p 



Y 



m 



Y 



where 



Y 



-k + l 

The discrete reverse Taylor remainder satisfies 
\eN-k{^^f{x + l),X,Y)\ 



(21) <max(2f~\l 



d.>0 JgN^ 
|J|<Af-fc 



\m 



{N-k-\3\+l) 
Y 



m 



{N-k-\J\+l),p^ 
Y 1 



Proof By definition ^ and (fT9l) . 

|5^_fc(AV(x + l),X,y)| <E[ 



< 



N,pE[ J2 ix+ji + ir]< 

o<ji<-<jN-k+i<y 



J] |A^+V(^ + l+ii 

o<ii<-<ijv-fc+i<i^ 

Y 

N,p 



LViV-A; + l 



{X + YY 



<max(2P-i,l)||/|U,p(E[XP]m 



y 



Y 
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where we have used in the last inequahty the estimations {X + < 2^ ^{X^ + Y^) if 
p > 1 and (X + Y)P < XP + YP Hp <1. Thus 1^ is proved. The inequahty 1^ follows 
from (El) and (EOl). □ 



Proposition 4.2 Let G N, p > and h G "Hp. Let Xj (i = 1, • • • , n) 6e a family of 
independent N-r.v.s with mean Aj > and up to {N + p + 1)*^ order moments; W = 
Xi + ■ ■ ■ + Xn- Let X* be an N-r.v. having Poisson Xi-zero biased distribution and 
independent ofW^^^ := W — Xi. Then the following estimations hold. 

1) WhenN = {), 
(22) 

n 

\e,{h)\ < m^x{2P~\im\\o,,Y.[E[{W^^^^)P]{E[Xf]+X^^ 

i=l 

2) When N > I, one has the recursive estimation: 
(23) 

n r 



eN{h)\<Y.Xi E mf;>{rnf)+m%^)\e^^\,\{^\'\fj,{x + l))\ 



i=l d>l JgN^ 
|J|<JV 

N 



fc=0 d>0 JeN^ 

|J|<Ar-fc 

+ max(2^-\ 1)||MU,, (E[(W»)^]m5^+^) + m%^'~^^^) 
Proof. We begin by the case when N = {). By 

n 

eo(/i) = 5] A, (E[A(iy W + XI + 1)] - E[A(VF« + X, + 1)]) 

i=l 
n 

= ^Wh{x + 1), W^'\X*) + eo(A(2; + 1), W^'^,Xi)) 

i=l 

n 

< max(2P-\l)||A||o,p^Ai{E[(t^«)P](m5;l +mg) + {m^'^ + m^^f)] 



i=l 



where the last inequality is by estimations (j20p and (j2ip . so (|22p follows. Combining in 
addition the recursive formula (llip . we obtain the inequality (I23p . □ 
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